Interchain monomer contact probability in two-dimensional polymer solutions 
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Using molecular dynamics simulation of a standard bead-spring model we investigate the den- 
sity crossover scaling of strictly two-dimensional (d = 2) self-avoiding polymer chains focusing on 
properties related to the contact exponent 82 set by the intrachain subchain size distribution. With 
R ~ N v being the size of chains of length N, the number njnt of interchain monomer contacts per 
monomer is found to scale as n lnt ~ l/N" 2 with v = 3/4 and 82 = 19/12 for dilute solutions and 
v = 1/d and 62 = 3/4 for TV 3> g(p) ~ 1/p 2 - Irrespective of the density p sufficiently long chains are 
thus found to consist of compact packings of blobs of fractal perimeter dimension d p = d — 82 = 5/4. 



As first suggested by de Gennes [l| , it is now generally 
accepted [2Hl3j that strictly two-dimensional (2D) self- 
avoiding polymer chains adopt compact and segregated 
conformations at high densities, i.e., the typical chain size 
R scales as 



R(N,p) ps {N/pY where v = 1/d = 1/2 



(1) 



with N being the chain length, p the monomer num- 
ber density, and d — 2 the spatial dimension. We as- 
sume here that monomer overlap and chain intersection 
are strictly forbidden Q. Compactness obviously does 
not imply Gaussian chain statistics nor does segregation 
of chains impose disklike shapes minimizing the contour 
perimeter of the chains. Focusing on dense 2D melts it 
has been shown in fact [H, fllT - [l3l ] that the irregular chain 
contours are characterized by a fractal perimeter of typ- 
ical length 



Nn-, r 



Rd P „ N d v /d with d d 



2 = 5/4 (2) 



where nj nt stands for the fraction of monomers of a chain 
interacting with monomers of other chains. The fractal 
line dimension d p is set by Duplantier's contact expo- 
nent 82 = 3/4 characterizing the size distribution of inner 
chain segments Q • Since the melt limit is experimentally 
difficult to realize for a strictly 2D layer [4|-|6| , the aim of 
this Letter is to show that eq holds more generally for 
all densities if the chains are sufficiently long. 

As described in Ref. [13] our numerical results are ob- 
tained by molecular dynamics simulations of monodis- 
perse, linear and highly flexible chains using the well- 
known Kremer-Grest bead-spring model 14|. Lennard- 
Jones (LJ) units are used throughout this paper and tem- 
perature is set to unity. Note that monomer overlap and 
chain crossings are made impossible by the excluded vol- 
ume interaction and the bonding potential; the chains 
are thus strictly 2D "self-avoiding walks" in the sense of 
Ref. 3. While our previous studies 11- 13| have focused 
on one melt density, p = 0.875, we scan now over a broad 
range of densities p as may be seen from FigureQ] Chains 
of length N = 16 up to N = 1024 monomers have been 
sampled for all p; for p — 0.125, p = 0.5 and p = 0.875 
we have even computed systems with N = 2048 fl5j ]. 
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FIG. 1: Average chain size R for self-avoiding polymers in 
d — 2 dimensions for a broad range of densities p and chain 
lengths N. The exponent v = vo = 3/4 expected for the dilute 
swollen chain limit is given by the dashed lines, the compact 
chain exponent v = 1/d by the bold lines. Main panel: Root- 
mean square end-to-end distance R c /£ and radius of gyration 
Rg/£ as a function of reduced chain length N/g where we 
set £ = l/p 3//2 and g = 1/p 2 . Inset: Replot focusing on the 
p-scaling for constant chain length N. 



That sufficiently long 2D polymer chains become in- 
deed compact for all densities, as stated by eqQ] is shown 
in Figure [1] presenting the overall chain size R as char- 
acterized by the root-mean square end-to-end distance 
R e and the radius of gyration R g [3]. As one expects 
[3, 0] , the typical chain size is found to increase with a 
power-law exponent v = i/q = 3/4 in the dilute limit 
(dashed lines) and with v — 1/d for larger chains and 
densities in agreement with various numerical [8j-|l0j and 
experimental studies 0, H, 0] ■ The main panel presents 
data for different densities as a function of chain length 
N which are successfully brought to collapse by reseat- 
ing the horizontal axis with the number of monomers 
g(p) ps p£ d ~ i/pV^orf-i) _ p-2 spanning the semidilute 
blob and the vertical axis by its size £ ~ g v ° ~ 1/p 3 / 2 . 
For simplicity, all prefactors related to the dilute limit 
are set to unity, e.g. g = 1/p 2 and £ = 1/p 3 / 2 [I3|- An 
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x = r/R(s) 



FIG. 2: Main panel: Subchain size distribution G(r, s) be- 
tween the monomers n — N/4+ 1 and m — n + s — 1 = 3-/V/4 
for N = 1024 and various densities using the same symbols 
as in the main panel of Figure Q] Motivated by the scaling 
in the dilute and melt limits, eq|3] we trace R d (s)G(r, s) vs. 
x — r/R(s) with R 2 {s) being the second moment of G(r,s). 
Inset: Interchain monomer pair distribution function giut(r) 
as a function of x = r/R s (N), revealing for i < 1 the same 
exponents 82 as the intrachain distribution G(r,s). 



equivalent data representation focusing on the chain size 
as a function of density taking as reference the dilute 
chain size Rq = N v ° is shown in the inset. We trace here 
R c /Ro as a function of the reduced density p/p* with 
p* = N/R d = 1/s/N being the crossover density. The 
typical chain size decreases indeed as R ~ with 
increasing density (bold line) as expected for compact 
chains, eqQ] Below we refer to N/g = (p/p*) 2 > 1 as 
the compact chain limit and to g — >■ as the melt limit. 

Obviously, the mere fact that v becomes 1/2 in the 
compact chain limit does not imply Gaussian chain statis- 
tics [3, This can be directly seen, e.g., from the sub- 
chain size distribution G(r, s) of the vector r = r_ m — r n 
between the monomers n and m = n + s — 1 on the 
same tagged chain presented in the main panel of Fig- 
ure [5] [18]. For the examples indicated we have chosen 
subchains of length s = N/2 in the middle of chains of 
length N — 1024 [l9j |. Only distances r > 2 are pre- 
sented to take off local physics related to the LJ beads. 
Using the measured second moment R 2 {s) of the distri- 
bution we trace R d (s)G(r, s) as a function of the reduced 
distance x = r/R(s). This is motivated by the scaling 

R d (s)G(r, s) = f(x) with f(x) ps x 02 for x < 1 (3) 

observed as expected in the dilute and melt limits where 
R(s) is the only relevant length scale. Note that f(x) is 
a cut-off function for x ^ 1 [l2j. The dashed and bold 
power-law slopes indicated in the figure correspond to the 
exponents 62 = 19/12 and 82 =3/4 predicted for both 
limits Q. (We remind that for Gaussian chains 62 = 0, 
i.e. G(r, s) decays monotonously with distance r.) Ob- 
viously, the scaling of G(r, s) becomes more intricate for 



semi-dilute densities where the blob size £ sets an ad- 
ditional length scale. For s >• g(p) 1 one expects 
to observe two power-law asymptotes, one scaling as 
G(r, s) ~ r 19 / 12 if the conformational properties within 
the blob are probed (r -C £) and one as G(r, s) ~ r 3 / 4 for 
5«r« R(s). Although much larger blobs are certainly 
warranted to demonstrate this unambiguously, this ex- 
pected behavior is qualitatively_ consistent with our data 
as may be seen for p — 0.125 [l9j. 

As shown in the inset of Figure [21 the 62 exponents 
do not only characterize intrachain properties but are 
also relevant for describing the interchain monomer con- 
tacts. We have plotted here the radial pair distribution 
function gi nt (r) between monomers from different chains 
which becomes unity (by normalization) for large reduced 
distances x — r/R g ^> 1. The scaling of the horizontal 
axis is again motivated by the density limits where the 
chain size R is the only relevant length scale. This scaling 
can be verified by comparing different chain lengths N at 
same density p (not shown). The power-law exponents 
8 2 for x 1 can be understood by generalizing a scaling 
argument given in Ref. Q stating that it is irrelevant of 
whether the neighborhood at r <C R around the reference 
monomer is penetrated by a long loop of the same chain 
or by another chain with a center of mass at a typical 
distance R. Hence I20B . 



7int 



(r)p fa G(r, s ps N) x N x R d p/N 



(4) 



where the second factor N stems from the fact that it is 
irrelevant which monomer of the second chain is probed. 
Note that the underlined term characterizing the prob- 
ability for having a second chain close to the reference 
chain drops out in the compact chain limit, eq[TJ Equa- 
tions [3] and [4] imply 



9bxt(r) 



l/N^ 2 for x < 1 



(5) 



as indicated by dashed and bold solid lines for, respec- 
tively, dilute and melt limits. 

In order to characterize more systematically the den- 
sity crossover we focus now on g- m t(r ps a) with a being 
one fixed monomeric distance for all densities. For the 
data presented in Figure [3] we have determined the typ- 
ical number m n t of interchain monomer contacts for a 
reference monomer by integrating g-mt(r)p up to a cut- 
off a = 1.56 2l|. Similar scaling results have been ob- 



tained by computing the fraction of monomers with at 
least one interchain contact within a distance a following 



Ref. 12| . Using a similar representation as in Figure [IJ 



we present the rescaled number of interchain contacts 
function of the reduced chain length N/g for 
different densities (main panel) and as a function of the 
reduced density p/p* for different chain lengths (inset). 
The rescaling used for the vertical axes is readily obtained 
by matching at N/g — (p/p*) 2 = 1 the power-law results 
implied by eq[S]for the dilute (v&2 — 19/16) and compact 
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FIG. 3: Average interchain monomer contact number ni nt vs. 
reduced chain length N/g (main panel) and reduced density 
p/p* (inset). The dashed lines correspond to the exponent 
19/16 for dilute chains, the full lines to u02 = 3/8 



vv 2 



for the compact limit. Confirming eqs [5] and [6] the observed 
successful data collapse is the central result of this work. 



(v92 =3/8) limits. This implies for the latter limit (solid 
bold lines) 



FIG. 4: Kratky representation of the intramolecular form fac- 
tor F(q). The Debye formula is indicated by the top thin 
line, the dilute swollen chain limit with F(q) ~ l/q 1 '" by 
the dashed line, the Porod scattering for a compact 2D ob- 
ject with smooth perimeter by the dash-dotted line and the 
power-law exponent —82 = —3/4 by the bold solid line. Our 
complete theoretical prediction for asymptotically long chains 
using the so-called Redner-des Cloizeaux approximation of 
G(r, s) allowing to fit the data starting from the Guinier 
regime for small wave vectors, is given by the thin solid line. 
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(6) 



where the first factor corresponds to the fraction of 
monomers of a given blob interacting with monomers 
from other chains given that the blob is at the chain 
perimeter. The second term stands for the fraction of 
blobs of a chain interacting with other chains and the suc- 
cessful scaling presented in Figure |31 thus demonstrates 
that in the compact chain limit the perimeter length L 
scales as L/£ w {N/g) 1 ' 92 /' 1 Since on the other 

hand L/£ « (i?/£) dp ~ (N/g) dp/d , this implies a fractal 
line dimension d p — d — 02=5/4, generalizing thus the 
numerical result obtained for the melt limit [12J to the 
broader compact chain limit. 

Although direct real-space visualizations of 2D poly- 
mer systems are experimentally feasible 0, 0, 0] , allow- 
ing thus a direct computation of G(r, s) and gintM, it is 
important to note that the exponent 62 = 3/4 for com- 
pact chains of blobs may in principle be obtained from 
an analysis of the intrachain molecular form factor F{q) 
with q being the wave vector 16[ . This point is addressed 
in Figured] where the rescaled form factor (F(q)/N)Q 2 is 
plotted as a function of the reduced wave vector Q = qR g 
for different densities p. Filled symbols refer to data for 
N = 2048, all other symbols to N = 1024. At variance 
to the indicated Debye formula for Gaussian chains a 
strong nonmonotonous behavior is revealed by our data, 
which approaches with increasing p and N a power-law 
exponent — 82 = —3/4 (bold solid line). This power law 
limit is implied by the fact that F(q) can be expressed 
as an integral over the Fourier transform G(q, s) of the 



intrachain size distribution G(r, s) [12|,ll8|]. Using general 
arguments which do not depend explicitly on the density 
or the blob size it has been shown [l2| that as a direct 
consequence of eq [3] the form factor for asymptotically 
long chains must scale as 



1 



1 



r 1 1 

for - < q < - (7) 



N ^ Q»2 Qd-d p 

where in the last step of the equation chain we have used 
dp = d — 82- We emphasize that this result is consis- 
tent with the well-known Porod scattering expected for 
gen eral compact objects of fractal surface dimension d p 
[1 31 ] . Obviously, much larger chains as computationally 
feasible at present are required to confirm numerically 
the predicted perimeter dimension using eq[7] 

In summary, we have investigated the density crossover 
scaling of strictly 2D self-avoiding polymer chains fo- 
cusing on properties related to the 62 contact exponent 
which characterizes not only (by definition) the intra- 
chain distribution G(r, s), but also the interchain pair dis- 
tribution function <?; nt (r) ~ r® 2 and the monomer contact 
number n^t ~ \/N" &2 . The established scaling for dense 



melts 0, [Tl]) [12| is found to remain valid for all densities 
if reformulated in terms of chains of blobs. As a conse- 
quence, the intramolecular form factor F(g) measuring 
the composition fluctuations at the fractal perimeters of 
the compact blob chains approaches with increasing den- 
sity and chain length a power-law asymptote implied by 
the generalized Porod scattering of compact objects of 



fractal perimeter dimension d p 



5/4. Note 



that our chains are sufficiently long to confirm d p using 



4 



the intermolecular binary contacts rii n t (Figure[3]) but not 
using the form factor. One reason for this is that eq[7| 
assumes asymptotically long chains where critical expo- 
nents other than # 2 may be ignored. The corrections 
for finite N/g need still to be worked out theoretically. 
However, since g- m t{r) or n lnt may be directly obtained 
experimentally by visualization of dilute and semidilute 
2D polymer solutions [1, [H, 0] > our numerical results sug- 
gest that the analysis of the latter properties may allow 
the experimental verification of the predicted exponents. 

Finally, we emphasize that at variance to most experi- 
mental systems 0,0,0 we have assumed in the presented 
work that the chains are flexible down to monomcric 
scales. As long as the persistence length remains smaller 
than the blob size, this should not alter the suggested 
scaling properties. In fact, since the blob size decreases 
very strongly with persistence length, i.e. N/ g increases, 
a finite rigidity might even speed up the convergence to 
the predicted asymptotic behavior. In order to bring our 
computational approach closer to experiment we are cur- 
rently computing systems with finite persistence length. 
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FIG. 5: TOC figure: Key scaling argument relating the in- 
trachain size distribution G(r, s ~ N) to the interchain pair 
distribution function gmt(r). It states that it is irrelevant of 
whether the neighborhood around the reference monomer n is 
penetrated by a long loop of the same chain as in panel (a) or 
by another chain with a center of mass at a typical distance 
R as shown in panel (b). 



